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Abstract 

We analyse the electric dipole moment of the neutron in the MSSM, induced by 
the renormalisation of the soft-susy breaking terms. We run the RGEs using two- 
loop expressions for gauge and Yukawa couplings and retaining family dependence. 
The fi and B parameters were determined by minimising the full one-loop Higgs 
potential, and we find that the neutron EDM lies in the range 10-33 < \dn\ < 
10-29 g cni. 
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Supersymmetric unified theories are the most promising candidates for physics be- 
yond the standard model in that they resolve the crucial gauge hierarchy problem of 
widely separated electroweak and grand unified scales, are the consequence of string the- 
ories and are favoured over non-supersymmetric unified theories by recent high precision 
measurements at LEP. In addition to the usual signatures of grand unification such as 
proton decay, neutrino masses, fermion mass relations and weak mixing angle prediction, 
supersymmetric unification is characterised by the resultant mass spectrum of the super- 
symmetric particles (squarks, sleptons, charginos, neutralinos) and the flavour- changing 
and CP-violating processes which arise as the renormalisation group equations (RGE) 
scale the physics from the unification scale Mjj ~ 10^^ GeV down to the electroweak 
scale. Of particular interest are the flavour changing neutral current transitions involving 
the quark-squark-gluino vertex, with their implications for rare B-decays and mixings 

and the non-removable CP- violating phases, with their implications for quark electric 
dipole moments (EDM) 0, ^ Q and for non-standard-model patterns of CP-violation 
in neutral B decays [Q, which result from the RGE scaling of the soft supersymmetry 
(SUSY) breaking scalar interactions in these models. 

CP-violation in the standard model (SM) arises from the single phase 6ckm in the 
Cabbibo-Kobayashi-Maskawa (CKM) mixing matrix relating the quark weak interaction 
and mass eigenstates and, in principle, this source of CP violation can accommodate the 
known CP properties of the kaon system. Non-zero quark (and lepton) electric dipole mo- 
ments are very sensitive probes of CP violation beyond the standard model [§] because, 
unlike the other observables of CP violation which are small because of the intergener- 
ational mixing angles of the CKM matrix, electric dipole moments are particularly sup- 
pressed by the chiral nature of the weak interaction and vanish at both one- and two-loop 
order in the Standard Model, resulting in quark EDMs of 

d'J-O{10-'') ecm. (1) 

At present the experimental bounds on quark EDMs are obtained indirectly from mea- 
surements of the neutron EDM. In the non-relativistic quark model the neutron EDM 
is 

4 1 

dn{quarks) = -dd - -du (2) 

and is of the same order as the u and (i-quark EDMs. However the neutron EDM is 
expected to be dominated by long distance effects such that 

d^^{LD) = 0(10-^2) e cm. (3) 

By contrast, non-zero supersymmetric phases, collectively denoted ^fj/f^, arising from 
the complexity of the soft SUSY breaking terms can generate quark EDMs at one-loop 
order, irrespective of generation mixing from diagrams involving gluino, chargino (or 
neutralino) exchange and mixing of right- and left-handed virtual squarks. However, for 
squark, gluino and chargino masses of order 100 GeV, these induced one-loop quark 
EDMs yield a neutron EDM which exceeds the experimental upper bound of 

rff^^ < 12 X 10-'^ e cm (4) 



unless these soft phases are constrained to satisfy 

■'soft 



cr < 0.01. (5) 
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The alternative scenario of not imposing any condition of smallness on these phases but 
instead making the supersymmetric scalar masses heavy enough ( of the order of 1 TeV) 
to suppress the EDMs has been considered by Kizukuri and Oshimo|0]. This scenario also 
has consequences for the relic density of the lightest supersymmetric particle 0. 

For supersymmetric unified models such as spontaneously broken = 1 supergravity 
with flat Kahler metrics ^ the resultant explicit soft SUSY breaking terms at the scale 
MsusY ~ Mjj of local SUSY breaking are quite simple and these generic phases 5fojf^ 
are reduced to just two phases Sj^^^ in addition to the usual Sckm- The most natural 
way of satisfying the experimental bound @ and ensuring that the squark, gluino and 
chargino masses are not much above the electroweak scale is to assume that these phases 
^A^B^^ vanish identically at the unification scale because of CP conservation in the SUSY 
breaking sector. Under these conditions the only explicit CP violation at the unification 
scale is in the flavour-dependent Yukawa coupling matrices which are required to have the 
structure necessary to reproduce the CKM mixing matrix at the electroweak scale under 
RGE scaling. However, as the RGEs for the soft SUSY breaking trilinear couplings Au^^^e 
and bilinear coupling B depend on the Yukawa couplings, inclusion of flavour mixing in the 
RGEs can lead to large RGE-induced CP- violating phases in the off-diagonal components 
of the couplings triggered, in particular, by the complexity of the large t-quark Yukawa 
coupling. 

The implications of such large phases for EDMs of quarks have been studied recently 
by Bertolini and Vissani |^ and Inui et al @] within a = 1 supergravity inspired 
minimally supersymmetric standard model (MSSM) in which the spontaneous breaking 
of the electroweak SU(2) x f/(l) symmetry is driven by radiative corrections. Bertolini 
and Vissani argue that the dominant induced EDM is that of the d-quark arising from 
the one-loop diagram involving chargino exchange and find 

dSUSY^^y^ ~ O{10'^^) e cm, (6) 

four orders of magnitude greater than the standard model prediction (|I|) but still satisfying 
the experimental upper bound (Q). Inui et a/ also find that the (i-quark EDM from chargino 
exchange is dominant but obtain the much larger value 

df^^{INUI) ~ 0(10-2^ - 10-29) e cm (7) 

which they ascribe to the inclusion of gaugino masses in their RGEs. 

Recently Dimopoulos and Hall have considered quark and lepton EDMs in a class of 
supersymmetric unified theories based on the gauge group SO(IO) where the unification of 
all quarks and leptons of a particular generation into a single 16 spinorial representation 
leads to non-removable CKM-like phases in the Yukawa couplings which, under RGE 
scahngs induced by a large t-quark Yukawa coupling, give rise to EDMs close to the 
experimental limits such that some regions of the parameter space of the minimal SO(IO) 
theory are excluded. 

The sensitivity of the quark EDMs to CP violation, and the fact that the EDMs of 
the Standard Model, the MSSM and GUT theories are nicely separated makes this an 
important window to physics at the unification scale. In addition the dipole moments in 
the MSSM are a minimum prediction of supersymmetry. Because of this, and also because 
of the numerical discrepancies between the calculations of Bertolini and Vissani and Inui 
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et al for the EDM of the (i-quark in the constrained MSSM, and the hmited nature of the 
free parameters chosen in both sets of calculations, we have undertaken a more detailed 
study of quark EDMs. 

We have used two loop evaluation of gauge and Yukawa couplings, rather than the 
one-loop RGEs as used for the existing EDM calculations, and minimised the full one- 
loop Higgs potential, including contributions from matter and gauge sectors. We do this 



following the very complete analyses of Kane et al |]T0[ and Barger et al , but retaining 



the full flavour dependence in the RGEs as we run them. We do not consider it necessary 



to describe the entire procedure since this is outlined in some detail in refs.[|TT], |T^, but 
we shall briefly consider some details of our analysis. 
The superpotential of the MSSM is given by 

W = KQIh^Ur + hdQ[H,DR + KL^H.En + fiH^eH^, (8) 

where generation indices are implied, and we define the VEVs of the Higgs fields (f i and 
^2) such that m„ = /i„f2, rrid = h^Vi and me = hf.Vi. Supersymmetry may be broken 
softly by generic mass-squared scalar terms, gaugino masses, and by 'trilinear' couplings 
of the form, 

5C = AuKqLh2UR + AdhdQihidR + A^hJliieR + Bfihieh2, (9) 

where again, generation indices are suppressed. In order to determine the Yukawa cou- 
plings at the weak scale, we first ran the Standard Model down to 1 GeV using the 
two- loop QED and three-loop QCD RGEs of Arason et al |jl2|. We then ran the full su- 
persymmetry RGEs up to the GUT scale (i.e. where the SU{2) and U{1) gauge couplings 
unified). For this we used two- loop RGEs for the gauge and Yukawa couplings flSl, and 



one-loop RGEs for everything else, in the DR scheme. Here we 'unify' by setting the 
strong coupling equal to the unified SU{2) and U{1) couplings. This neglects the effects 
of thresholds at the GUT scale which depend on the precise details of the GUT theory, 
and tends to give values of (~ 0.126 for a top mass of 174 GeV) at the weak scale 
which are a little on the high side |T^. At the GUT scale there are three parameters 



which we set by hand, the common scalar mass mg, the common gaugino mass mi/2, and 
the common trilinear coupling A. Some of this degeneracy (for example of the gaugino 
masses) is motivated by the presumed existence of a GUT theory, and some by minimal 
supergravity (together with the assumption that the effects of renormalisation between 
the GUT and Planck scales is small). From the point of view of determining the effects 
of renormalisation of the pure MSSM on the electric dipole moments, degeneracy of these 
parameters is the natural assumption. We then run the entire theory back down to the 
weak scale. We determined the mass eigenvalues and eigenstates at the relevant physical 
scale, Q = m{Q). At the same time we retained the full supersymmetric spectrum in 
the running theory, and only decoupled states in the running of the gauge and Yukawa 
couplings. Given the value of tan/5 and the sign of /i, it is possible to minimise the ef- 
fective potential using the tadpole equations of ref. |]ll| in the running theory. This is a 
valid procedure if we wish to minimise the one-loop effective potential only, and avoids 
the need to match Lagrangians at each particle threshold ||T0|, The minimisation was 
done at Q = rutop (which was taken to be 174 GeV throughout), leaving the canonical, 
hybrid, four- dimensional, parameter space (mo, m.1/2. A, tan/3) in addition to sign(/i). 
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However, since all the contributions from matter and gauge sectors were included in the 
minimisation, the vacuum expectation values of the two neutral Higgs fields (or equiva- 
lently the values obtained for B and /i) should be insensitive to the momentum scale at 
which they are evaluated |T^. This was indeed found to be the case. The whole process 



was then iterated a number of times. Generally, the procedure converges very rapidly 
(within a few iterations), and we accurately recover the entire spectrum given in ref. [|rT[ 
for a wide range of parameters^ 

The mass matrices (in the super-KM basis) were diagonalised numerically to yield the 
required mass eigenstates and diagonalisation matrices as follows, 

squarks : V^M?Vq = m? 
neutralinos : V^MatVat = m^o 

charginos : Uf^McVc = ""^x*' i^^) 



and the following constraints applied. 

In addition to experimental constraints (we adopt those used in ref. [pi) 
18 GeV, m^± > 47 GeV, niho > 44 GeV, m/,± > 44 GeV, tua > 21 GeV, m-g > 141 GeV, 
rriy > 43 GeV, rriq > 45 GeV) we insisted that the minimum was stable in the sense that 
the Higgs and squark mass-squareds were positive. We also required that the minimum 
which we obtained was global, and that there were no other minima which may have 
broken colour or charge. The constraints 

Al < 3hl{ml + ml + ml) 
Al < 3hl{ml + ml + ml) 

A? < 3hUml + ml + ml) (11) 

provide a coarse indication of this [0. We also insisted that the lightest supersymmetric 
partner was the neutralino, and finally we required that the process converged (i.e. that 
our choice of parameters was not too close to any fixed points). 

The diagonalisation matrices appear in trilinear couplings between the quarks and the 
heavy supersymmetric scalar bosons and fermions, in particular squarks and charginos 
or gluinos. It is CP violation (i.e. non-zero phases) in these matrices, at the interaction 
vertices of the diagrams shown in figs(la,lb,2), which may induce a non-zero EDM. As 
discussed in ref.[^, we find that such phases are indeed induced into the A-terms by the 
running of the RGEs, and hence into the diagonalisation matrices. 

Having established this fact, let us consider how the EDM arises. For completeness, 
we wish to include the w-quark contribution (which is usually neglected), and so we shall 
briefly re-examine the EDM calculation. In doing so we also hope to gain a little insight 
into the CP violating nature of the EDM. First focus on a gaugeless Lagrangian with two 
fermionic fields, one scalar field, and a single cubic coupling, 

C = ipiii'j''df,-mi)ipi+i!2{h''d^-m2)i'2 
+ 19^01' -mj 101' 

+ («V^2L^li? + blp2R^lL) 0* + (a>i^?/'2L + &*V'iL^2r) 0- (12) 



We would like to thank P. Ohmann for discussions 
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The ijji field is the light quark whose EDM we would like to calculate, the scalar field, 0, 
represents the squark fields, and ip2 is the heavy fermion field (not another quark). For 
any particular quark, we may choose that basis in which the mass parameters are real. 
A CP transformation on this Lagrangian shows that, in order to have CP violation, the 
phases of a and b must be different (a common phase may be absorbed into the definition 
of 4>). Now consider the self energy graphs in fig(la,b). In addition to giving mass and 
wave function renormalisation, these diagrams also induce non-local terms which may be 
obtained by performing a derivative expansion 

AC = D[d^i^,j,d^ij,L-mli^,jiij,L)G{x) 

+ D* (d^i^.r.d'^ip.R - m^V^i^Viii) G{x) + ... (13) 

where the dots represent terms which are higher order in momentum, and where 

G{x) = ^3^(l-a;2 + 2xlogx) 
D = a*b/ (3271^7712) 

X = ml/ ml- (14) 

With CP violation D is complex. The EDM appears when we now introduce electromag- 
netic interactions whilst keeping this expansion gauge invariant by introducing covariant 
derivatives, 

^^^0 {df, + iq^A^)(f) (15) 

with qi = q2 + q^- When D is complex, one can anticipate an electric dipole moment from 
A£, of 

d^qiG{x)lm{D). (16) 

In fact, when the heavy fermion is electrically neutral (like the gluino of fig(la)), this is 
the one-loop contribution to the dipole moment. In general there is an additional gauge 
invariant contribution to the Lagrangian, coming from the heavy fermion charge. To 
determine this we must resort to the usual one-loop diagram shown in fig(2). It is found 
to be of the form 

A£' = -q2F^,H{x) (Re(D)V^ia'^>i - ilm{D)'^^a>'''-f''il;,) , (17) 

where, 

Ffxi/ = d^Ay — di^A^, 

^[7^71, (18) 



and where 



2 

= — - X + xXogx) . (19) 

(1 — Xf 
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The total EDM of a quark coming from chargino/squark loops is then, 

d = Im(L>) [qiG{x) - q2H{x)] . 
For each of the quarks this gives, 



(20) 



dd = -lm{Dd)Fd{x) 
2e 

du = -yIm(L'„)F„(a;), 



(21) 



where we have defined the functions, 

1 



(1-X)3 
1 

(1 - xf 



5 - 12a; + + 2a;(2 - 3a;) logo; 



2 — 6x + Ax^ + x{l — 2>x) \ogx 



(22) 



The above analysis generalises in a straightforward manner. The coupling constants a and 
h become matrices a-ij and fejj, with i,j running over the appropriate mass eigenstates. 
For the chargino contributions we find 



dd = 



3 327r2 ^ 



V 



^ 2 e 

d„ = -- 



3 327r2 



VuFd 



2 

m ± 

X, 



2 

TTlj 
g 

2 

xf 



K^hd 



RL 
T 
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(23) 



RL 
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where we are using the down-quark diagonal basis, and where K is the CKM matrix. For 
the gluino contributions we find. 



dd 
du 



-Im 



LR/ 11 



2eas 
dTrrrig 



■Im 



VuG 



m:- 



2lK 



(24) 



Li?/ 11 



In order to present our results, we choose points generated at random in the parameter 
space given by 



< 


mo 


< 1 TeV 


< 


mi/2 


< 1 TeV 


-1 < 


A 


< 1 TeV 


< 


tan/3 


< 20. 



(25) 

In practice, values higher than these seldom satisfy all the criteria detailed above (i.e. 
they imply fine-tuning). The region below the low tan/3 fixed point is excluded. For a 
top-quark of mass 174 GeV, A — {) GeV and mi/2 — rriQ — 150 GeV, this was found to 
be at tan/3 = 2.2. 
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The modulus of the neutron EDM is plotted against tan/? in fig(3). There is a slight 
tendency for it to be positive, and the largest values occur for negative /i and positive A. 
Clearly the value of tan f3 dominates the EDM of the neutron, and we see the approx- 
imately linear behaviour for large tan/3 coming from the increased down-quark Yukawa 
couplings. Phases feed into the off-diagonal elements of A^, especially into A^is. The 
EDM becomes smaller as we approach the low tan/5 fixed point, since the top Yukawa 
coupling dominates the running. Here the gluino contribution to the up-quark can be the 
dominant contribution. Elsewhere however, the down-quark, chargino diagram is nearly 
always dominant. No obvious pattern emerges with the other three parameters. The 
value of the neutron EDM in the MSSM is much less than the value of 10~^^ indicated 
in 0], because values of the A parameter as large as those used in ref.0], give problems 
with colour or charge breaking minima, or do not lead to a solution for /i and B on 
minimisation. We find that the expected range for the EDM is therefore 

10"^^ < < 10-29ecm. (26) 

It is expected that future developments will push the experimental bound on the 
neutron EDM down from eq.(|D to O(10~^^)e cm but, as the present calculations indicate, 
will still not provide a test of the constrained MSSMQ. A more remote possibility is the 
direct measurement of the t-quark EDM using either tt decay correlations in e~^e~ — >■ tt 
\[% or tt production via photon-photon fusion using linearly polarised photons generated 



by Compton back-scattering of laser light on electron or positron beams of linear e~^e 



or e~e~ colliders In this analysis, we found that the t-quark EDM is usually larger 
than the neutron EDM by a factor of 3-5, a slight improvement but still unlikely to be 
measured in the forseeable future, at least in the MSSM. 



^We thank K. Green for discussions 
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Figure Captions 

Figure 1 : Quark self energy diagrams involving (a) gluino and (b) chargino exchange. 

Figure 2 : SUSY contribution to quark EDM from chargino-photon coupling. 

Figure 3 : The modulus of Neutron EDM (in units of 10"'^^ e cm) as a function of tan j3 for 
random choices of mo, mi/2 and \A\'m the range (0, 1) TeV. 
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